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A model of nuclear fusion consisting of a wave packet impinging into a well located between square 
one dimensional barriers is treated analytically. The wave function inside the well is calculated 
exactly for the assisted tunneling induced by a perturbation mimicking a constant electric field with 
arbitrary time dependence. Conditions are found for the enhancement of fusion. 
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I. INTRODUCTION 

Quantum tunneling explains the transmission of par- 
ticles or clusters of particles through regions that are 
energetically forbidden. A time-honored success of the 
quantum tunneling model is the explanation of a de- 
cay lifetimes of unstable nuclei by Gurney and Condon 
and Gamowp], 0]. More recent revisions of the tunneling 
paradigm in the nuclear decay context, have reaffirmed 
its validity. [1 

In previous works we investigated the tunneling of a 
one dimensional metastable state between barriers, ex- 
cited by a time dependent potential. [j| , 

In ref.Q we found analytical expressions for the as- 
sisted tunneling processes. The acceleration of the decay 
of a metastable state was found to be determined by the 
poles of the unperturbed wave functions in the complex 
energy plane. 

In the present work we generalize the approach of 
and find exact solutions to the assisted tunneling of a 
wave packet impinging from a region outside the range 
of the potentials, into a square well lying between square 
barriers. This setup provides a simplified model for the 
fusion of nuclei. 

Fusion time scales are of the order of the transit time 
of the nuclei past each other and differ from a decay 
lifetimes markedly. Nuclear decay lifetimes are typically 
very long in comparison to natural nuclear times. On this 
basis, it may be argued that fusion can not be affected in 
the same manner as found for the case of a decay. Nev- 
ertheless, in the presence of a long range perturbation, 
the tunneling process is sensible to longer time scales. 
Consequently, an external agent can affect the inter pen- 
etration the nuclei. 

However, for long times, it is not possible to follow the 
decay process numerically. Q . Hence, analytical formulas 
are vital for nuclear fusion as well as for decay. This is 
the motivation for the present effort. 

In the next section we review the model and extend it 
the case of fusion. In section 3 we apply the results to 
the penetration of a wave packet into a region between 



square barriers. By inspecting the formulas we are able 
identify the relevant parameter space for the enhance- 
ment of fusion. 



II. ANALYSIS OF ASSISTED TUNNELING 

The Schrodinger equation for a one dimensional system 
consisting of a square well between square barriers is [l5[ 
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The inclusion of a square well between the barriers is 
needed to model nuclear fusion. 

The even and odd stationary states of eq. |T]) for ener- 
gies below the barrier strength 7 are 



Xe,o( x ) 
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where 
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cos(qx); \x\ < Xo 
ip e {k) = { A lG K W + Bie~ K W ■ d > \x\ > x (3) 
Cicos(kx) + sign(x)Disin(kx); \x\ > d 

where k and q are defined in eq. (|26p . and sign(x) = 1,-1 
for x > 0,x < 0. 



sin(qx); \x\ < xq 
<p {k) = { sign{x){A 2 e^ x \ + B 2 e~^ x \); d > \x\ > x (4) 
sign(x)C2Cos(kx) + D2sin(kx); \x\ > d 

, and we have extracted a factor of y/n from the normal- 
izations for convenience. The labels e, o refer to the even 
or odd character of the wave functions 

The set of even-odd functions is orthonormal and 
complete. [ill [I2J The all important normalization fac- 
tors of eq.((2]), and the amplitudes of eqs. (|3I4|) . are shown 
in appendix A. 
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We subject the system to a time dependent spatially 
linear perturbation, higher powers in space can be dealt 
similarly 



For an initial (t = 0), wave packet representing the in- 
cident nucleus carrying mean momentum qq and centered 
at ro, ^ q „,r {x,t — 0), the solution of eq. (fTT|) reads 



V(x,t)=nx G(t) 



(5) 



with [i a coupling constant. For the case of a spatially 
constant time-harmonic electric field of intensity Eq in- 
teracting with a nucleus of charge Z |e| we have 



[i G(t) = Z \e\ Eq sin(u>t) 
Applying a unitary transformation 
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C = J G(t) dt 



(7) 



The Schrodinger equationfl]) including the perturbation 
of eq.© becomes 



l ^~t = Y^dtf +7(e(a:o + d-N)e(H- a: o))* 



v e(x -\x\)$ + i((t) 



/J 



<9$ 

dx 



(8) 



Eq. ijSJ) is solved by expanding the wave function in the 
complete set of even and odd states of the unperturbed 
Schrodinger equation © 



/•OO 



e 2 



dk (9) 



Taking advantage of the superposition integral evaluated 
in the appendix A of adapted now to the case of a 
potential well between the barriers 



dx 



n e (k)n (k) 



S(k-k') 



(10) 



with q defined in eq. (|26[) of Appendix A, the time evolu- 
tion of the amplitudes c ej0 of eq.([9]) becomes 



c (k) 
c e {k) 



i« q C 



to n e (k) n a (k) 
M Q C 



c e (k) 



to n e (k) n {k) 
a dot denoting a time derivative. 



c a (k) 



(11) 



C e (k,t) 

c (k,t) 



Y{t) 



A e cos(z) + A sin(z) 
A cos(z) — A e sin(z) 
fiqY 



to n e n a 



C(t') dt' 



(12) 



where 



A e 
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n e 

Do 
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^q ( ur Q Xo dx 



(13) 



III. 



ASSISTED TUNNELING INTO A WELL 
BETWEEN BARRIERS 



The most important contribution to eq.© in the re- 
gion of the well is due to the zeros of the normalization 
factors n e of eq.© displayed in eqs. (|24l25[) of appendix 
A.1,1 

Figure 1 shows the inverse of the normalizations of 
eq.© for to = 15000Afe^,x = 7/m, d = 12/to, 7 = 
10M eV, Vq = 40MeV, corresponding to an le O nucleus 
for energies below the barrier. The spikes in figure 1 
are due to the extreme closeness of the minima of the 
normalization factors to their complex zeros. 

For eigenenergies above the barrier, the inverse of the 
amplitudes show a smooth spectrum of peaks of order 
one in height and width. These peaks will contribute a 
negligible amount at long times due to the strong oscil- 
lations of the time dependent exponential factors of the 
wave functions, and are henceforth omitted. 
The number of poles is directly related to the width 
and height of the barrier. Very few minima and conse- 
quently very few poles show up below the barrier energy 
for lighter nuclei, as compared to the a decay process 
from heavier nuclei for which the repulsive barriers are 
wider and taller. [6j 

The normalization factors around their minima can be 
cast in the form|| 



k 



Aj,(e,o) {k Z ~ fc 3 ,( e ,o) 2 ) 2 + Pj,(e,o) (14) 



where j enumerates the zeros of n e ,n . 

The poles appear in pairs located symmetrically above 
and below the real momentum axis[6(. The stationary 
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FIG. 1: Normalizations of eq.Q around the minima at k 
ko, with Ak = k — fco, seetext. 



wave functions of eqs. (|3l4[) consist of incoming and out- 
going parts that have poles in the upper and lower halves 
of the complex energy plane respectively [Hj]. 

The minima of the normalization factors are extremely 
sharp. It was necessary to apply a multiple precision 
numerical package to locate them Q with a maximum of 
50 decimal digits. The imaginary parts of the poles are 
orders of magnitude smaller than the real parts. (l6j 

Eq.@ can be performed now by contour integration 
in the complex energy plane, closing the contour from 
below to insure convergence. The negative imaginary 
part of each pole induces a decaying exponential when 
the dominant time dependence is due to the phase of 
the unperturbed wave function. The perturbation can 
induce both decaying and growing behaviors. In 
we were interested in the acceleration of the decay. Here, 
we seek an enhancement of the rate of penetration into 
the well region, for which a different parameter space will 
be needed. 

After performing the complex energy plane contour in- 
tegration, the wave function in the region between the 
barriers becomes 



$ = $mt + ®ext 
&ext = £l,e + £2,0 + £3,0 + £,i.e 



(15) 



The splitting into an external and internal contribution 
is necessary because the outer part of the wave function 
vanishes at a pole. The treatment of the contour inte- 
gration has to be modified accordingly. 

The rather involved expressions of the wave functions 
of eqs. fTS"]) are spelled out in appendix B. 

Using the lar ge a rgument expansion of the hypergeo- 
metric functions in eqs. (129H33p and eqs. (|34j|37p of Ap- 
pendix B, the long time t — ► oo behavior of eq. (fT5|) in the 
region of the well |x| < Xq becomes 



^ Vt 

3 

Hkf k] Y(t) kf k) t 
m \q \ m 



(16) 



where we have used the connection between the even and 
odd amplitudes at the poles @ 



no{k od ) 2 0t d 



= k 
= k 



2 

e,j 
2 



(17) 



kj, kj denote the real and imaginary parts of the poles, 

I qo I = y2 m Vq, and Cj(x) is a prefactor that depends 
on the specific features of the initial wave packet. 

Inspection of eq. (|16p reveals that the perturbation can 
either accelerate the tunneling process of slow it down 
depending on the sign of F(t). This is also true for the 
assisted tunneling corresponding to a decay. In [j| , Q we 
focused on assisted decay for which F(i) was demanded 
to be negative. 

It follows from eq. (QJJ|) that the perturbation will assist 
the fusion process when 



(xY(t) 



> t 



m\ 



(18) 



independently of the pole structure. 

To fulfill eq. (fT5)l . Y(i) has to be positive definite. 
Using eqs. $fn§ ,G(t) = sin(tt t), with fl ^ 0, 



Y(t) 



(Q t - sin(Q t)) 

«2 



(19) 



Y(t) > is satisfied for O t » 1. 
For long times, eq. lfl8|) becomes 



n < 



JJ_ 



(20) 



It is now possible now to estimate the parameters 
needed for assisted fusion. 

Consider low density and temperature totally ionized 
oxygen nuclei approximated as an ideal gas. For a tem- 
perature of 100°^ the velocity of the oxygen nuclei is 
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v « 400 . The inter nuclei distance at P = O.lPa 

sec 

is w 2400A The average time between collisions is 
St « 6 10" 10 sec. If we take ft £t = 10, eq.® becomes 

M > 10 M (21) 

Inserting the fusion parameters used in figure 1 into 
eq.([2"Tj) we find 

/x > 600^^ (22) 
cm 

or equivalently an electric field amplitude in eq.© of 

around \E \ « 8 10 7 -^. 

cm 

This electric field appears quite large but not unreach- 
able, especially in light of the fact that the corresponding 
frequency of eg. ([20]) is fi w YIGhz and the perturbation 



can be applied in ultrashort pulses of nanosecond dura- 
tion. 

For angular frequencies smaller than the value pre- 
scribed by eq. ([20|) . the tunneling into the well of a packet 
located far away outside the well, will increase as com- 
pared to the unassisted case. 

Undoubtedly, the actual fusion problem is more com- 
plicated than the simplified model of a packet impinging 
on a well between barriers, especially because of nuclear 
structure aspects. However, the analysis of tunneling 
provided here suggests that the enhancement of fusion 
by means of external time dependent agents is possible. 
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UNPERTURBED WAVE FUNCTIONS 



2 e\ e 2 qk 
+ e 2 K 2 si s 2 — e\ q k C\ s 2 

The amplitudes of the unperturbed wave functions of + e 2 q k ci c 2 + e 2 k k si c 2 

eqs. (|3l4l) are + e\ n k s\ c 2 — e\ q k c\ c 2 ) 
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Do 



1 



I 2 ,22 

(e 2 q k c\C2 + e 2 K si C2 

2 



2 ei e2 g 
— e\n 2 Si c 2 + e\ q n C\ c 2 
+ e\ q k C\ s 2 + e 2 k k si s 2 
+ e\ k k S\ s 2 — e\ q k ci s 2 



(25) 



ft' 

q 

ei 

e-2 

si 

Cl 

si 



= 



to 7 - 



k 2 



y/2 to (k 2 + V ) 



cos{k d) 
sin(k d) 
cos(q Xq) 
sin(q xq) 



(26) 



APPENDIX B: WAVE FUNCTION IN THE REGION 
BETWEEN THE BARRIERS 



The even and odd amplitudes of eq. (|T5|) are written as 



Ann — 



A e '° 4- A e '° 



A e '° - - 



/OO 
Q(|x| - d)^ q0iro Xe,o dx 
-OG 

Q(d- \x\)^g 0tro Xe,o dx 



(27) 



Splitting of the contributions of the inner and outer re- 
gions is needed because the outer part of the wave func- 
tions vanishes at the poles. The finite part of the time 
dependent wave function for the outer regions has to be 
evaluated by a slightly different method than the contri- 
bution of the inner part. 

For the sake of exemplification, we consider an initial 
Gaussian wave packet carrying mean momentum qo and 
centered at ro, 



% Ot r o (x,t = 0) = 

<i> = iq ( x -r )- {x ' X 9 ° )2 (28) 



Inserting eq. ([27|) into the wave function of eqs. lfTS]) for 
the region between the barriers, the contour integration 
around the poles yield 



E ( n i)2 

j n=0 v ' 



Q< 



V^F A™' coafe, e x) e v 



2 y // 2f3j e \j 
oo 

E E« 



{-xloY 



'{n\) 2 (2n + l) (n + 1) 



V^A oC0S (fc,, o i)e"> 2 F(t)V 2 



where 



X j,(e,o) = M 9j,(e,o) 



2 m n (kj^ efi ))y/Pj,(e,o) 



~ k j,(e,o)(- 



,A 2 



2~m 



Qj,(e,o) = \J \ k j,(e,o)) 2 + 2 to V 



oo /_ y2 \n 

E E ' ■ '' 



(n!) 2 (2 n+1) 



j n=0 

o cos(kj e x') e 



2 to n 2 (kj. e ) y/2/3j :e \ jie 



E E 



(-*£„)" 

(n!) 2 (2 n + 1) 



j n=0 

H Y(t) q 3 . ^K nt ) e 



2 to n 2 (fc 3 ,o) >/2/3j i0 Aj )0 



E E< 



n=0 



(n!) 2 (2n + l)(n+ 1) 



-M 2 ^M 2 g 2 e V^A™*szn(% : e a;)e t ' 



■=■3,0 



4 v / 2/3 J - e A J - e TO 2 n4(fc iie 

OO 

E E QooL 7n+7 



(-*L) n 



j n=0 

Qoo = VttAq sin(k jt0 x) e l 



2 ^/2/3j Aj i0 



(-*L) n 



Q oe = — M q Jje n/ttA™' sin{kj, e x) 



(29) 



(30) 



(31) 



(32) 
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2 to n 2 {k^ e ) y/2/3~X~ 



(n!) 2 (2 n+1) 



j n=0 

-/i qj M \pKK e sin(kj^ x) e 1 
■ 

''i.fi 



2 to nl(k j:0 ) y/2f3 j:0 X j:0 



(33) 



•s3,o 



^ oo 



\\3 



j n=0 v 

yfirA e xt sin{kj }0 x) e v 
2 



(36) 



^ eo 



ex I 



2-, 2-, / 2 n3 

VttA^ 4 cos(kj. e x) e v 



(2n 



E 



t (4 n)! 

(2 n!) 2 (2 n+1) 

^ q jt0y /nAl xt cos{kj, x) e v 
2 to n 2 (fc i:0 ) ^A^o 



i n=0 



(34) 



(35) 



(-^, e ) n (4 n)! 
(2 n!) 2 (2 n+1) 



j n=0 

/i y(t) qj^eV^ex^e sin(kj, e x) e v 
2 to nl(kj^) J\j, e 



(37) 



where fcj. e is the complex momentum at the pole number j 
of the even unperturbed wave function, Y(t) is defined in 
eq. (112|l , and A and (3 correspond to the expansion around 
a pole of ea.(|14p. The wave functions of eqs. (|29H3"3")) and 
eqs. (|34H37j) , can be expressed in terms of standard hyper- 
geometric, Struve, and Bessel functions, fbj] 



